Impact of public health policy and mobility change on transmission potential of severe acute respiratory syndrome coronavirus 2 in Rhode Island, March 2020 - November 2021.

To study the SARS-CoV-2 transmission potential in Rhode Island (RI) and its association with policy changes and mobility changes, the time-varying reproduction number, Rt, was estimated. The daily incident case counts (16 March 2020, through 30 November 2021) were bootstrapped within a 15-day sliding window and multiplied by Poisson-distributed multipliers (λ = 4, sensitivity analysis: 11) to generate 1000 estimated infection counts, to which EpiEstim was applied to generate Rt time series. The median Rt percentage change when policies changed was estimated. The time lag correlations were assessed between the 7-day moving average of the relative changes in Google mobility data in the first 90 days, and Rt and estimated infection count, respectively. There were three major pandemic waves in RI in 2020-2021: spring 2020, winter 2020-2021 and fall-winter 2021. The median Rt fluctuated within the range of 0.5-2 from April 2020 to November 2021. Mask mandate (18 April 2020) was associated with a decrease in Rt (-25.99%, 95% CrI: -37.42%, -14.30%). Termination of mask mandates on 6 July 2021 was associated with an increase in Rt (36.74%, 95% CrI: 27.20%, 49.13%). Positive correlations were found between changes in grocery and pharmacy, Rt retail and recreation, transit, and workplace visits, for both Rt and estimated infection count, respectively. Negative correlations were found between changes in residential area visits for both Rt and estimated infection count, respectively. Public health policies enacted in RI were associated with changes in the pandemic trajectory. This ecological study provides further evidence of how non-pharmaceutical interventions and vaccination slowed COVID-19 transmission in RI.

Excess mortality in Ukraine during the course of COVID-19 pandemic in 2020-2021.

In this work, the COVID-19 pandemic burden in Ukraine is investigated retrospectively using the excess mortality measures during 2020-2021. In particular, the epidemic impact on the Ukrainian population is studied via the standardized both all-cause and cause-specific mortality scores before and during the epidemic. The excess mortality counts during the pandemic were predicted based on historic data using parametric and nonparametric modeling and then compared with the actual reported counts to quantify the excess. The corresponding standardized mortality P-score metrics were also compared with the neighboring countries. In summary, there were three "waves" of excess all-cause mortality in Ukraine in December 2020, April 2021 and November 2021 with excess of 32%, 43% and 83% above the expected mortality. Each new "wave" of the all-cause mortality was higher than the previous one and the mortality "peaks" corresponded in time to three "waves" of lab-confirmed COVID-19 mortality. The lab-confirmed COVID-19 mortality constituted 9% to 24% of the all-cause mortality during those three peak months. Overall, the mortality trends in Ukraine over time were similar to neighboring countries where vaccination coverage was similar to that in Ukraine. For cause-specific mortality, the excess observed was due to pneumonia as well as circulatory system disease categories that peaked at the same times as the all-cause and lab-confirmed COVID-19 mortality, which was expected. The pneumonias as well as circulatory system disease categories constituted the majority of all cases during those peak times. The seasonality in mortality due to the infectious and parasitic disease category became less pronounced during the pandemic. While the reported numbers were always relatively low, alcohol-related mortality also declined during the pandemic.

SARS-CoV-2 Transmission in Alberta, British Columbia, and Ontario, Canada, December 25, 2019, to December 1, 2020.

OBJECTIVE: This study aimed to investigate coronavirus disease (COVID-19) epidemiology in Alberta, British Columbia, and Ontario, Canada. METHODS: Using data through December 1, 2020, we estimated time-varying reproduction number, Rt, using EpiEstim package in R, and calculated incidence rate ratios (IRR) across the 3 provinces. RESULTS: In Ontario, 76% (92 745/121 745) of cases were in Toronto, Peel, York, Ottawa, and Durham; in Alberta, 82% (49 878/61 169) in Calgary and Edmonton; in British Columbia, 90% (31 142/34 699) in Fraser and Vancouver Coastal. Across 3 provinces, Rt dropped to ≤ 1 after April. In Ontario, Rt would remain < 1 in April if congregate-setting-associated cases were excluded. Over summer, Rt maintained < 1 in Ontario, ~1 in British Columbia, and ~1 in Alberta, except early July when Rt was > 1. In all 3 provinces, Rt was > 1, reflecting surges in case count from September through November. Compared with British Columbia (684.2 cases per 100 000), Alberta (IRR = 2.0; 1399.3 cases per 100 000) and Ontario (IRR = 1.2; 835.8 cases per 100 000) had a higher cumulative case count per 100 000 population. CONCLUSIONS: Alberta and Ontario had a higher incidence rate than British Columbia, but Rt trajectories were similar across all 3 provinces.

SARS-CoV-2 Transmission Potential and Policy Changes in South Carolina, February 2020 - January 2021.

INTRODUCTION: We aimed to examine how public health policies influenced the dynamics of coronavirus disease 2019 (COVID-19) time-varying reproductive number (R t ) in South Carolina from February 26, 2020, to January 1, 2021. METHODS: COVID-19 case series (March 6, 2020, to January 10, 2021) were shifted by 9 d to approximate the infection date. We analyzed the effects of state and county policies on R t using EpiEstim. We performed linear regression to evaluate if per-capita cumulative case count varies across counties with different population size. RESULTS: R t shifted from 2-3 in March to <1 during April and May. R t rose over the summer and stayed between 1.4 and 0.7. The introduction of statewide mask mandates was associated with a decline in R t (-15.3%; 95% CrI, -13.6%, -16.8%), and school re-opening, an increase by 12.3% (95% CrI, 10.1%, 14.4%). Less densely populated counties had higher attack rates (P < 0.0001). CONCLUSIONS: The R t dynamics over time indicated that public health interventions substantially slowed COVID-19 transmission in South Carolina, while their relaxation may have promoted further transmission. Policies encouraging people to stay home, such as closing nonessential businesses, were associated with R t reduction, while policies that encouraged more movement, such as re-opening schools, were associated with R t increase.

Age-Stratified Model to Assess Health Outcomes of COVID-19 Vaccination Strategies, Ghana.

We assessed the effect of various COVID-19 vaccination strategies on health outcomes in Ghana by using an age-stratified compartmental model. We stratified the population into 3 age groups: <25 years, 25-64 years, and ≥65 years. We explored 5 vaccination optimization scenarios using 2 contact matrices, assuming that 1 million persons could be vaccinated in either 3 or 6 months. We assessed these vaccine optimization strategies for the initial strain, followed by a sensitivity analysis for the Delta variant. We found that vaccinating persons <25 years of age was associated with the lowest cumulative infections for the main matrix, for both the initial strain and the Delta variant. Prioritizing the elderly (≥65 years of age) was associated with the lowest cumulative deaths for both strains in all scenarios. The consensus between the findings of both contact matrices depended on the vaccine rollout period and the objective of the vaccination program.

Comparative epidemiology of five waves of COVID-19 in Mexico, March 2020-August 2022.

BACKGROUND: The Mexican Institute of Social Security (IMSS) is the largest health care provider in Mexico, covering about 48% of the Mexican population. In this report, we describe the epidemiological patterns related to confirmed cases, hospitalizations, intubations, and in-hospital mortality due to COVID-19 and associated factors, during five epidemic waves recorded in the IMSS surveillance system. METHODS: We analyzed COVID-19 laboratory-confirmed cases from the Online Epidemiological Surveillance System (SINOLAVE) from March 29th, 2020, to August 27th, 2022. We constructed weekly epidemic curves describing temporal patterns of confirmed cases and hospitalizations by age, gender, and wave. We also estimated hospitalization, intubation, and hospital case fatality rates. The mean days of in-hospital stay and hospital admission delay were calculated across five pandemic waves. Logistic regression models were employed to assess the association between demographic factors, comorbidities, wave, and vaccination and the risk of severe disease and in-hospital death. RESULTS: A total of 3,396,375 laboratory-confirmed COVID-19 cases were recorded across the five waves. The introduction of rapid antigen testing at the end of 2020 increased detection and modified epidemiological estimates. Overall, 11% (95% CI 10.9, 11.1) of confirmed cases were hospitalized, 20.6% (95% CI 20.5, 20.7) of the hospitalized cases were intubated, and the hospital case fatality rate was 45.1% (95% CI 44.9, 45.3). The mean in-hospital stay was 9.11 days, and patients were admitted on average 5.07 days after symptoms onset. The most recent waves dominated by the Omicron variant had the highest incidence. Hospitalization, intubation, and mean hospitalization days decreased during subsequent waves. The in-hospital case fatality rate fluctuated across waves, reaching its highest value during the second wave in winter 2020. A notable decrease in hospitalization was observed primarily among individuals ≥ 60 years. The risk of severe disease and death was positively associated with comorbidities, age, and male gender; and declined with later waves and vaccination status. CONCLUSION: During the five pandemic waves, we observed an increase in the number of cases and a reduction in severity metrics. During the first three waves, the high in-hospital fatality rate was associated with hospitalization practices for critical patients with comorbidities.

Late surges in COVID-19 cases and varying transmission potential partially due to public health policy changes in 5 Western states, March 10, 2020-January 10, 2021.

OBJECTIVE: This study investigates the SARS-CoV-2 transmission potential in North Dakota, South Dakota, Montana, Wyoming, and Idaho from March 2020 through January 2021. METHODS: Time-varying reproduction numbers, R t , of a 7-day-sliding-window and of non-overlapping-windows between policy changes were estimated utilizing the instantaneous reproduction number method. Linear regression was performed to evaluate if per-capita cumulative case-count varied across counties with different population size or density. RESULTS: The median 7-day-sliding-window R t estimates across the studied region varied between 1 and 1.25 during September through November 2020. Between November 13 and 18, R t was reduced by 14.71% (95% credible interval, CrI, [14.41%, 14.99%]) in North Dakota following a mask mandate; Idaho saw a 1.93% (95% CrI [1.87%, 1.99%]) reduction and Montana saw a 9.63% (95% CrI [9.26%, 9.98%]) reduction following the tightening of restrictions. High-population and high-density counties had higher per-capita cumulative case-count in North Dakota on June 30, August 31, October 31, and December 31, 2020. In Idaho, North Dakota, South Dakota and Wyoming, there were positive correlations between population size and per-capita weekly incident case-count, adjusted for calendar time and social vulnerability index variables. CONCLUSIONS: R t decreased after mask mandate during the region's case-count spike suggested reduction in SARS-CoV-2 transmission.

An ensemble n-sub-epidemic modeling framework for short-term forecasting epidemic trajectories: Application to the COVID-19 pandemic in the USA.

We analyze an ensemble of n-sub-epidemic modeling for forecasting the trajectory of epidemics and pandemics. These ensemble modeling approaches, and models that integrate sub-epidemics to capture complex temporal dynamics, have demonstrated powerful forecasting capability. This modeling framework can characterize complex epidemic patterns, including plateaus, epidemic resurgences, and epidemic waves characterized by multiple peaks of different sizes. We systematically assess their calibration and short-term forecasting performance in short-term forecasts for the COVID-19 pandemic in the USA from late April 2020 to late February 2022. We compare their performance with two commonly used statistical ARIMA models. The best fit sub-epidemic model and three ensemble models constructed using the top-ranking sub-epidemic models consistently outperformed the ARIMA models in terms of the weighted interval score (WIS) and the coverage of the 95% prediction interval across the 10-, 20-, and 30-day short-term forecasts. In our 30-day forecasts, the average WIS ranged from 377.6 to 421.3 for the sub-epidemic models, whereas it ranged from 439.29 to 767.05 for the ARIMA models. Across 98 short-term forecasts, the ensemble model incorporating the top four ranking sub-epidemic models (Ensemble(4)) outperformed the (log) ARIMA model 66.3% of the time, and the ARIMA model, 69.4% of the time in 30-day ahead forecasts in terms of the WIS. Ensemble(4) consistently yielded the best performance in terms of the metrics that account for the uncertainty of the predictions. This framework can be readily applied to investigate the spread of epidemics and pandemics beyond COVID-19, as well as other dynamic growth processes found in nature and society that would benefit from short-term predictions.

Investigating COVID-19 transmission and mortality differences between indigenous and non-indigenous populations in Mexico.

OBJECTIVES: Indigenous populations have been disproportionately affected during pandemics. We investigated COVID-19 mortality estimates among indigenous and non-indigenous populations at national and sub-national levels in Mexico. METHODS: We obtained data from the Ministry of Health, Mexico, on 2,173,036 laboratory-confirmed RT-PCR positive COVID-19 cases and 238,803 deaths. We estimated mortality per 1000 person-weeks, mortality rate ratio (RR) among indigenous vs. non-indigenous groups, and hazard ratio (HR) for COVID-19 deaths across four waves of the pandemic, from February 2020 to March 2022. We also assessed differences in the reproduction number (Rt). RESULTS: The mortality rate among indigenous populations of Mexico was 68% higher than that of non-indigenous groups. Out of 32 federal entities, 23 exhibited higher mortality rates among indigenous groups (P < 0.05 in 13 entities). The fourth wave showed the highest RR (2.40). The crude HR was 1.67 (95% CI: 1.62, 1.72), which decreased to 1.08 (95% CI: 1.04, 1.11) after controlling for other covariates. During the intense fourth wave, the Rt among the two groups was comparable. CONCLUSION: Indigenous status is a significant risk factor for COVID-19 mortality in Mexico. Our findings may reflect disparities in non-pharmaceutical (e.g., handwashing and using facemasks), and COVID-19 vaccination interventions among indigenous and non-indigenous populations in Mexico.

Parameter identifiability and optimal control of an SARS-CoV-2 model early in the pandemic.

We fit an SARS-CoV-2 model to US data of COVID-19 cases and deaths. We conclude that the model is not structurally identifiable. We make the model identifiable by prefixing some of the parameters from external information. Practical identifiability of the model through Monte Carlo simulations reveals that two of the parameters may not be practically identifiable. With thus identified parameters, we set up an optimal control problem with social distancing and isolation as control variables. We investigate two scenarios: the controls are applied for the entire duration and the controls are applied only for the period of time. Our results show that if the controls are applied early in the epidemic, the reduction in the infected classes is at least an order of magnitude higher compared to when controls are applied with 2-week delay. Further, removing the controls before the pandemic ends leads to rebound of the infected classes.

Transmission Dynamics of COVID-19 in Ghana and the Impact of Public Health Interventions.

This study characterized COVID-19 transmission in Ghana in 2020 and 2021 by estimating the time-varying reproduction number (Rt) and exploring its association with various public health interventions at the national and regional levels. Ghana experienced four pandemic waves, with epidemic peaks in July 2020 and January, August, and December 2021. The epidemic peak was the highest nationwide in December 2021 with Rt ≥ 2. Throughout 2020 and 2021, per-capita cumulative case count by region increased with population size. Mobility data suggested a negative correlation between Rt and staying home during the first 90 days of the pandemic. The relaxation of movement restrictions and religious gatherings was not associated with increased Rt in the regions with fewer case burdens. Rt decreased from > 1 when schools reopened in January 2021 to < 1 after vaccination rollout in March 2021. Findings indicated most public health interventions were associated with Rt reduction at the national and regional levels.

SARS-CoV-2 transmission potential and rural-urban disease burden disparities across Alabama, Louisiana, and Mississippi, March 2020 - May 2021.

PURPOSE: To quantify and compare SARS-CoV-2 transmission potential across Alabama, Louisiana, and Mississippi and selected counties. METHODS: To determine the time-varying reproduction number Rt of SARS-CoV-2, we applied the R package EpiEstim to the time series of daily incidence of confirmed cases (mid-March 2020 - May 17, 2021) shifted backward by 9 days. Median Rt percentage change when policies changed was determined. Linear regression was performed between log10-transformed cumulative incidence and log10-transformed population size at four time points. RESULTS: Stay-at-home orders, face mask mandates, and vaccinations were associated with the most significant reductions in SARS-CoV-2 transmission in the three southern states. Rt across the three states decreased significantly by ≥20% following stay-at-home orders. We observed varying degrees of reductions in Rt across states following other policies. Rural Alabama counties experienced higher per capita cumulative cases relative to urban ones as of June 17 and October 17, 2020. Meanwhile, Louisiana and Mississippi saw the disproportionate impact of SARS-CoV-2 in rural counties compared to urban ones throughout the study period. CONCLUSION: State and county policies had an impact on local pandemic trajectories. The rural-urban disparities in case burden call for evidence-based approaches in tailoring health promotion interventions and vaccination campaigns to rural residents.

An investigation of spatial-temporal patterns and predictions of the coronavirus 2019 pandemic in Colombia, 2020-2021.

Colombia announced the first case of severe acute respiratory syndrome coronavirus 2 on March 6, 2020. Since then, the country has reported a total of 5,002,387 cases and 127,258 deaths as of October 31, 2021. The aggressive transmission dynamics of SARS-CoV-2 motivate an investigation of COVID-19 at the national and regional levels in Colombia. We utilize the case incidence and mortality data to estimate the transmission potential and generate short-term forecasts of the COVID-19 pandemic to inform the public health policies using previously validated mathematical models. The analysis is augmented by the examination of geographic heterogeneity of COVID-19 at the departmental level along with the investigation of mobility and social media trends. Overall, the national and regional reproduction numbers show sustained disease transmission during the early phase of the pandemic, exhibiting sub-exponential growth dynamics. Whereas the most recent estimates of reproduction number indicate disease containment, with Rt<1.0 as of October 31, 2021. On the forecasting front, the sub-epidemic model performs best at capturing the 30-day ahead COVID-19 trajectory compared to the Richards and generalized logistic growth model. Nevertheless, the spatial variability in the incidence rate patterns across different departments can be grouped into four distinct clusters. As the case incidence surged in July 2020, an increase in mobility patterns was also observed. On the contrary, a spike in the number of tweets indicating the stay-at-home orders was observed in November 2020 when the case incidence had already plateaued, indicating the pandemic fatigue in the country.

The doubling time analysis for modified infectious disease Richards model with applications to COVID-19 pandemic.

In the absence of reliable information about transmission mechanisms for emerging infectious diseases, simple phenomenological models could provide a starting point to assess the potential outcomes of unfolding public health emergencies, particularly when the epidemiological characteristics of the disease are poorly understood or subject to substantial uncertainty. In this study, we employ the modified Richards model to analyze the growth of an epidemic in terms of 1) the number of times cumulative cases double until the epidemic peaks and 2) the rate at which the intervals between consecutive doubling times increase during the early ascending stage of the outbreak. Our theoretical analysis of doubling times is combined with rigorous numerical simulations and uncertainty quantification using synthetic and real data for COVID-19 pandemic. The doubling-time approach allows to employ early epidemic data to differentiate between the most dangerous threats, which double in size many times over the intervals that are nearly invariant, and the least transmissible diseases, which double in size only a few times with doubling periods rapidly growing.

Time Reallocation and the Cost and Benefit of School Closures during an Epidemic

School closures are an important public health intervention during epidemics. Yet, the existing estimates of policy costs and benefits overlook the impact of human behavior and labor market conditions. We use an integrated assessment framework to quantify the public health benefits and the economic costs of school closures based on activity patterns derived from the American Time-Use Survey (ATUS) for a pandemic like COVID-19. We develop a policy decision framework based on marginal benefits and costs to estimate the optimal school closure duration. The results suggest that the optimal school closure depends on how people reallocate their time when schools are closed. Widespread social distancing behavior implemented early and for a long duration can delay the epidemic for years, buying time for the development of pharmaceutical interventions and yielding substantial net benefits. Conversely, school closure, with behavior targeted to adjust only to the school closure, is unlikely to provide substantial delay or sufficient net benefits to justify closing schools for pathogen control.

Transmission Dynamics and Short-Term Forecasts of COVID-19: Nepal 2020/2021.

Nepal was hard hit by a second wave of COVID-19 from April-May 2021. We investigated the transmission dynamics of COVID-19 at the national and provincial levels by using data on laboratory-confirmed RT-PCR positive cases from the official national situation reports. We performed 8 week-to-week sequential forecasts of 10-days and 20-days at national level using three dynamic phenomenological growth models from 5 March 2021-22 May 2021. We also estimated effective and instantaneous reproduction numbers at national and provincial levels using established methods and evaluated the mobility trends using Google's mobility data. Our forecast estimates indicated a declining trend of COVID-19 cases in Nepal as of June 2021. Sub-epidemic and Richards models provided reasonable short-term projections of COVID-19 cases based on standard performance metrics. There was a linear pattern in the trajectory of COVID-19 incidence during the first wave (deceleration of growth parameter (p) = 0.41-0.43, reproduction number (Rt) at 1.1 (95% CI: 1.1, 1.2)), and a sub-exponential growth pattern in the second wave (p = 0.61 (95% CI: 0.58, 0.64)) and Rt at 1.3 (95% CI: 1.3, 1.3)). Across provinces, Rt ranged from 1.2 to 1.5 during the early growth phase of the second wave. The instantaneous Rt fluctuated around 1.0 since January 2021 indicating well sustained transmission. The peak in mobility across different areas coincided with an increasing incidence trend of COVID-19. In conclusion, we found that the sub-epidemic and Richards models yielded reasonable short-terms projections of the COVID-19 trajectory in Nepal, which are useful for healthcare utilization planning.

Geospatial Variability in Excess Death Rates during the COVID-19 Pandemic in Mexico: Examining Socio Demographic, Climate and Population Health Characteristics.

OBJECTIVES: This study examined how socio-demographic, climate and population health characteristics shaped the geospatial variability in excess mortality patterns during the COVID-19 pandemic in Mexico. METHODS: We used Serfling regression models to estimate all-cause excess mortality rates for all 32 Mexican states. The association between socio-demographic, climate, health indicators and excess mortality rates were determined using multiple linear regression analyses. Functional data analysis characterized clusters of states with distinct excess mortality growth rate curves. RESULTS: The overall all-cause excess deaths rate during the COVID-19 pandemic in Mexico until April 10, 2021 was estimated at 39.66 per 10 000 population. The lowest excess death rates were observed in southeastern states including Chiapas (12.72) and Oaxaca (13.42), whereas Mexico City had the highest rate (106.17), followed by Tlaxcala (51.99). We found a positive association of excess mortality rates with aging index, marginalization index, and average household size (P < 0.001) in the final adjusted model (Model R2=77%). We identified four distinct clusters with qualitatively similar excess mortality curves. CONCLUSION: Central states exhibited the highest excess mortality rates, whereas the distribution of aging index, marginalization index, and average household size explained the variability in excess mortality rates across Mexico.

Investigating the first stage of the COVID-19 pandemic in Ukraine using epidemiological and genomic data.

The novel coronavirus SARS-CoV-2 was first detected in China in December 2019 and has rapidly spread around the globe. The World Health Organization declared COVID-19 a pandemic in March 2020 just three months after the introduction of the virus. Individual nations have implemented and enforced a variety of social distancing interventions to slow the virus spread, that had different degrees of success. Understanding the role of non-pharmaceutical interventions (NPIs) on COVID-19 transmission in different settings is highly important. While most such studies have focused on China, neighboring Asian counties, Western Europe, and North America, there is a scarcity of studies for Eastern Europe. The aim of this epidemiological study is to fill this gap by analyzing the characteristics of the first months of the epidemic in Ukraine using agent-based modelling and phylodynamics. Specifically, first we studied the dynamics of COVID-19 incidence and mortality and explored the impact of epidemic NPIs. Our stochastic model suggests, that even a small delay of weeks could have increased the number of cases by up to 50%, with the potential to overwhelm hospital systems. Second, the genomic data analysis suggests that there have been multiple introductions of SARS-CoV-2 into Ukraine during the early stages of the epidemic. Our findings support the conclusion that the implemented travel restrictions may have had limited impact on the epidemic spread. Third, the basic reproduction number for the epidemic that has been estimated independently from case counts data and from genomic data suggest sustained intra-country transmissions.

A Large-Scale COVID-19 Twitter Chatter Dataset for Open Scientific Research-An International Collaboration.

As the COVID-19 pandemic continues to spread worldwide, an unprecedented amount of open data is being generated for medical, genetics, and epidemiological research. The unparalleled rate at which many research groups around the world are releasing data and publications on the ongoing pandemic is allowing other scientists to learn from local experiences and data generated on the front lines of the COVID-19 pandemic. However, there is a need to integrate additional data sources that map and measure the role of social dynamics of such a unique worldwide event in biomedical, biological, and epidemiological analyses. For this purpose, we present a large-scale curated dataset of over 1.12 billion tweets, growing daily, related to COVID-19 chatter generated from 1 January 2020 to 27 June 2021 at the time of writing. This data source provides a freely available additional data source for researchers worldwide to conduct a wide and diverse number of research projects, such as epidemiological analyses, emotional and mental responses to social distancing measures, the identification of sources of misinformation, stratified measurement of sentiment towards the pandemic in near real time, among many others.

Harnessing testing strategies and public health measures to avert COVID-19 outbreaks during ocean cruises.

To ensure the safe operation of schools, workplaces, nursing homes, and other businesses during COVID-19 pandemic there is an urgent need to develop cost-effective public health strategies. Here we focus on the cruise industry which was hit early by the COVID-19 pandemic, with more than 40 cruise ships reporting COVID-19 infections. We apply mathematical modeling to assess the impact of testing strategies together with social distancing protocols on the spread of the novel coronavirus during ocean cruises using an individual-level stochastic model of the transmission dynamics of COVID-19. We model the contact network, the potential importation of cases arising during shore excursions, the temporal course of infectivity at the individual level, the effects of social distancing strategies, different testing scenarios characterized by the test's sensitivity profile, and testing frequency. Our findings indicate that PCR testing at embarkation and daily testing of all individuals aboard, together with increased social distancing and other public health measures, should allow for rapid detection and isolation of COVID-19 infections and dramatically reducing the probability of onboard COVID-19 community spread. In contrast, relying only on PCR testing at embarkation would not be sufficient to avert outbreaks, even when implementing substantial levels of social distancing measures.